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We study the validity of the Born-Oppenheimer approximation in chaotic dynamics. Using nu- 
merical solutions of autonomous Fermi accelerators, we show that the general adiabatic conditions 
can be interpreted as the narrowness of the chaotic region in phase space. 
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I. INTRODUCTION 

The Born-Oppenheimer approximation is a method 
which deals with a coupled system of heavy and light 
objects [ij. From atomic physics to nuclear physics, it 
has played an important role. Moreover, the importance 
of this method gets greater as the necessity of the adi- 
abatic control of quantum states arises in the quantum 
information science 

Nevertheless, the criteria of the validity of the method 
is not well defined yet, even still controversial in some 
aspects 0, Therefore, more theoretical work is needed 
through applications to specific models. 

In this paper, we focused on the effect of chaotic 
dynamics in the Born-Oppenheimer approximation 
method. Since the Born-Oppenheimer approximation 
basically deals with coupled systems, there is always 
the possibility of non-integrability generating chaotic 
dynamics 0. For this study, we choose the Fermi ac- 
celerator as a model of Hamiltonian chaos [ij. 

Originally, the Fermi accelerator has been proposed 
to explain the origin of fast cosmic rays and developed 
further by Ulam[7|. However, this model has had more 
importance as a standard Hamiltonian system to exhibit 
chaotic dynamics. Because of its simple composition and 
yet rich dynamics, it has been intensely studied both 
classically^ 0] and quantum mechanically 0, [13, flit . 
The feature of harmonically driven barriers is frequently 
encountered in various physical situations [13, [i3|j hence 
it provides very useful physical insights in spite of its 
simple structure. Furthermore, physical features of the 
model have been realized experimentally in an atomic 
system [3]. 

Besides the considerations of various regimes, several 
varieties of the Fermi accelerator have also been pro- 
posed, including relativistic models, of the type which 
embodies thermodynamic considerationsjlSj and me- 
chanical Fermi accelerator with damping [l6|. Our model 
can however be contrasted with the previous ones in 
terms of energy conservation. While the general Fermi 
accelerator has a periodic external driving force, the two 
moving objects in our system simply exchange energy 
with each other, but always maintaining the total en- 
ergy conserved. ( In the relativistic regime, the system 
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FIG. 1: Schemetical diagram of the autonomous Fermi accel- 
erator 



can have a constant energy effectively iii spite of the driv- 
ing, but it is conditionally given. See [l5|) Accordingly, 
this model may be called an 'Autonomous Fermi Accel- 
erator' (AFA), and the property of energy conservation 
enables us to calculate the exact eigenstates. With this 
set of eigenstates, we analyze the dynamics of the AFA 
and derive the condition for the validity of the Born- 
Oppenheimer approximation in the chaotic regime. 

This paper is organized as following. In Sec. II, we 
describe the autonomous Fermi accelerator schematically 
and present the analysis of its classical dynamics by phase 
space portraits. The phase space of the system shows un- 
der what conditions the chaotic dynamics sets in. Then 
we apply the Born-Oppenheimer approximation to the 
system in Sec. III. We then compare the approximated 
calculation with the exact numerical solution which is 
obtained in Sec. III. Through this comparison, the con- 
nection between the conventional criteria of the Born- 
Oppenheimer approximation and the chaotic dynamics 
can be derived. 



II. CLASSICAL DYNAMICS 

The autonomous Fermi accelerator consists of a free 
particle confined between two infinite potential barriers. 
Both barriers are not penetrable and one of them has a 
finite mass and moves in a harmonic potential, while the 
other one is fixed. The schematical description is shown 
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FIG. 2: Phase space portraits of the autonomous Fermi ac- 
celerator. The plots are constructed by recording dynamical 
phases at each collision between the light particle and the 
heavy wall : (a) , (b) Plots of dynamical phases of a light par- 
ticle, (c), (d) Same plots of a center of masses. The mass 
ratio(77) and the total energy(_B) for each plot are assigned as 
(r; = 0.01, E = 12.0) for (a) and (c), and (ry = 0.1, E = 12.0) 
for (b) and (d). 



As Fig. [2] shows, the phase space of the system is quite 
clearly divided into two parts, a stochastic and a regular 
part. The border between them is conditioned by the 
existence of the double collision. Here, the double colli- 
sion means the situation that the light particle is hit by 
the heavy wall twice until it reaches the fixed wall. The 
presence of this type of collision diverges the light par- 
ticle's trajectory drastically, and results in the stochas- 
tic dynamics. Apparently, it can be a necessary condi- 
tion for the double collision that the heavy wall is faster 
than the light particle. Fig. [5] shows such tendency well. 
The stochastic region is located around the equator of 
the sphere in the phase space, while the regular island is 
around the pole. 

We also mention a structural detail in phase space: 
Around the border between the chaotic and the regular 
region there exists "Stochastic Web Structure" [17)] . How- 
ever, this structure is ignorable in the regime of wave 
mechanical parameters in this study. 

The structural properties of the phase space, including 
the portion of the regular and chaotic region, is controlled 
mainly by the mass ratio(77). By the comparison of (a) 
and (b) in Fig. [2l we can confirm the tendency that the 
regular island is shrinked as r/ is increased. Also the phase 
space is affected by the total energy (i<^), but the portion 
of regular island is not changed as much as by the mass 
ratio, according to our numerical study. 
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The Hamiltonian of the autonomous Fermi accelerator 
is given by 



H{P,p,R,r) 
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where M, P and R are a mass, momentum and postion 
variable for the oscillating heavy wall, while m, p and r 
are those for the free particle. We set the fixed wall at 
R = and the minimum point of the harmonic potential 
at i?o = 5. Also, we assume that every collision between 
the particle and the wall is elastic. 

First, we analyze the classical dynamics of this system 
by Poincare surface of sections. The surface of section 
can be constructed stroboscopically by sampling the dy- 
namical status of either the particle or the wall at each 
collision. The obtained results with different parameters 
are presented in Fig. [2l The right column of the figure 
shows the phase space plots of the light particle and the 
left one shows those of the center of masses. 

In this figure, we can notice that the feature of energy 
conservation is vividly shown; namely that all dynamical 
phases are restricted on a spherical surface given by 
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III. BORN-OPPENHEIMER APPROXIMATION 

If the system is quantized, the total Hamiltonian be- 
comes 
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and the Schrodinger equation and the boundary con- 
dition are given by 



F*(r, R) = E<^{r, R) 



*(r = 0, i?) = < r, r) = 
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Application of the Born-Oppenheimer approximation 
means to fix all parameters which are given instanta- 
neously by the heavy object. In our system, this can 
be done by fixing the distance between the walls, to be 
i?, at every moment. Then, a set of basis states can be 
constructed as follows: 



^{R,r)=Y,UR)4>n{r,R) 
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where 
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Using the completeness of the basis set above, the 
Schrodinger equation ^ can be reduced in the foUow- 
ing way: 
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where e = {hmrY /2mR^ . 
If Eq. ([7]) is inserted into 



|), one is able to expand the 



equation with the mass-ratio, rj{= m/M) as a parameter. 
By taking up to the first order of this expansion, the 
following equation is obtained, 
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Consequently, the governing equation is reduced to a 
simple 1-dimensional equation of motion with a station- 
ary bounding potential and the third term of the equation 
can be thought of as a pressure generated by the frequent 
collisions of the light particle. 

According to the theory of Hamiltonian chaos, there is 
no possibility of chaos if we consider the dynamics at the 
pure classical level. However, if we take the level transi- 
tion into account, the effective potential can be thought 
of to be time-dependent, and it seems more reasonable 
as the counterpart of the classical chaotic system. Along 
this line of reasoning, if an eigenmode is corresponding 
to chaotic dynamics, we can make an assumption that 
the lowest order of approximation would fail to describe 
the eigenmode. The validity of this assumption, will be 
discussed in the next section through a discussion of the 
exact solution which is numerically obtained. 



IV. NUMERICAL STUDY 

As mentioned in the introduction, the exact solutions 
of the system are also available. To calculate them, the 
following transform is performed. 

R — X , r = rj y . (10) 
Then the Schrodinger equation becomes 



{-^V^ + \Kx-xof)^{x,y) = E^{x,y) . (11) 
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FIG. 3: Configurational Plots of Eigenmodes at r; = 0.01. : 

(a) Regular mode : E = 12.38 (c) Chaotic mode : E = 12.29. 

(b) , (d) : Corresponding Husimi distributions to (a) and (c), 
respectively 
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FIG. 4: Configurational Plots of Eigenmodes at 77 = 0.1: (a) 
Regular mode : E = 12.22 (c) Chaotic mode : E = 12.19. 
(b), (d) : Corresponding Husimi distributions to (a) and (c), 
respectively. 



As the result of the transformation, the 1-dimensional 
dynamics of the system becomes a 2-dimensional dynam- 
ics in a triangular well, the three sides of which are two 
infinite potential barriers along the x-axis and y — tan^a;, 
and one harmonic potential barrier to the x-direction. 
Here, the angle 9 between the two infinite barriers is given 
by 

d = tan^i ^ , (12) 
where 77 is the mass ratio(= m/M). It should be noted 
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that 9 would reach 7r/2 as rj approaches infinity. In this 
work, we are concerned only with the normal situation 
of Born-Oppenheimer approximation, namely that M is 
much larger than m. 

Then the boundary condition is given as 



= 0) = , 
'^{x,y = iSiuOx) = , 
^{x ^ oo,y) ^ . 



(13) 



Using a numerical approach, we obtain all possible 
eigenvalues and corresponding eigenmodes with two dif- 
ferent mass-ratio parameters, r] = 0.01 and 0.1. As 
the numerical algorithm, 'Finite Element Method' 18] is 
used. 

Fig.[3]and Fig.|4]show configurational plots and Husinii 
distributions of eigenmodes at 77 = 0.01 and ry = 0.1, re- 
spectively. Using the normal derivative of configurational 
mode distributions along y — r]x, we can calculate the 
Husimi distributions, and these are corresponding to the 
Poincare surface of section of a center of mass. Hence, 
the corresponding dynamics for each mode can be con- 
firmed by comparison of the Husimi distributions with 
the classical phase spaces in Fig. [5J In Fig. [3] and IH (a) 
and (b) of each figure shows a mode which is localized 
the regular regime in the phase space, and (c) and (d) 
shows ones in the chaotic regime. 

From the obtained computational results, the validity 
of the Born-Oppenheimer approximation is investigated. 
For this purpose, the Fourier analysis is performed to find 
the eigenmode composition of the light particle. 
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The result of the analysis shows the different features 
of distributions, depending on the underlying dynamics. 
For 7] = 0.01, the distribution is well localized regardless 
to the dynamics (see Fig. [5l^a).). For 77 = 0.1, however, 
the distribution shows the interesting features depending 
on the dynamics (Fig.[5Ub)). The mode which is placed in 
the chaotic region shows a broad distribution, with tails 
reaching the regular region, but the one in the regular 
region still has a well localized distribution. 

Considering the scheme of the Born-Oppenheimer ap- 
proximation, modes with the high localization are ex- 
pected to be well-approximated. To confirm this, we 
extract the wave function of the moving wall with the 
highest level (n) in each distribution in the following way. 
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Fig. [5] presents a comparison between the wave func- 
tion extracted from the exact wave function in Fig. [3Jc) 
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FIG. 5: (a) The analysis of the modes in Fig. E] (r? = 0.01) 
under the instantaneous basis of the light particle, (b) The 
same one for Fig. [l] In each figure, the upper distribution is 
for modes in the chaotic region and the lower one is for ones 
in the regular region, (c), (d) The extracted wave functions of 
the moving waIl(solid line) and the approximated wave func- 
tions(dashed line) 



and the one approximated by the Eq. ([9]). At 77 = 0.01, 
the extracted wave function agrees very well with the ap- 
proximation, regardless to where in phase space region a 
mode is situated (Fig. [D^c)). In contrast, the mode in 
the chaotic region for ry = 0.1 does not show a good 
agreement with the approximation, while the one in the 
regular region still abides by the aforementioned approx- 
imation (Fig. [Sjd)). Of course, the disagreement is due 
to the level transitions during the dynamics. 

Now, it is necessary to consider the general criteria for 
the variation of Born-Oppenheimer approximation [2^ 
2ll,[2|, which is given as follows. 



h\ < n\m > \^\En- E„ 



(16) 



This inequality simply guarantees that the heavy ob- 
ject is moving slowly enough to prevent level transitions. 
If we apply this inequality to our system, the condition 
of this inequality is not changed much by the dynamical 
status. In other words, the ratio between the left term 
and the right term in Eq. (|16p is given as around 0.01 for 
77 = 0.01 and around 0.2 for 77 = 0.1, no matter whether 
a mode is chaotic or regular. 

Considering the near-intcgrability in the case of regular 
dynamics, it is obvious that the simple approximation 
works well in the regular region, because we can define 
good quantum numbers for each degree of freedom. 

In addition, we can find the dynamical interpretation 
of the inequality (fTB|) . In our system, the satisfaction of 
guarantees that the chaotic region is small enough 
in comparison to Planck's constant h. Therefore, it is 
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impossible that more than two modes are situated in the 
chaotic region. Accordingly, no transition is possible. 



V. CONCLUSION 

In this work, we study the validity of the Born- 
Oppenheimer approximation in the presence of chaotic 
dynamics. Especially, the criteria under which the ap- 
proximation is working is interpreted from the dynami- 
cal point of view. We apply the approximation method 
to the autonomous Fermi accelerator which has a clear 
division of regular dynamics and chaotic dynamics in 
its phase space. In the regular region, the approxima- 
tion is obviously well suited to calculate the eigenmodes, 
whereas it is just conditionally suited in the chaotic re- 
gion. However, when the dynamics of the system satisfies 
the general condition for the approximation, that is, the 
evolution of the system is slow enough to prevent the 
level transitions, the area of the chaotic region becomes 
small enough in comparison to Planck's constant. Thus, 
the chaotic dynamics can not be resolved by wave func- 
tions. Thereby, the Born-Oppenheimer approximation 
can approximate the wave functions very well in spite of 
the chaotic nature of the dynamics. 

It is important to remark that in the region where 
the Born-Oppenheimer approximation is working, we 
can apply another approximation to our system. Then 
we can follow the procedure of Ref. to construct 
the density matrix and the reduced density matrix of 
the interesting subsystem (light particle or oscillating 
wall). We leave this for a future work. We also remark 
that it would be quite instructive to study higher order 
corrections to the Born-Oppenheimer approximation 
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leave this for a future study. 
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